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We investigate the problem of metric fluctuations in the presence of the vacuum fluctuations 
of matter fields and critically assess the usual assertion that vacuum energy implies a Planckian 
cosmological constant. To this end, a new stochastic classical approach to the quantum fluctuations 
of spacetime is developed. The work extends conceptually Boyer's random electrodynamics to a 
theory of random gravity but has a considerably richer structure for inheriting nonlinearity from 
general relativity. Attention is drawn to subtleties in choosing boundary conditions for metric 
fluctuations in relation to their dynamical consequences. We point out that those compatible with 
the observed Lorentz invariance must allow for spontaneous conformal fluctuations, in addition to 
stochastic gravitational waves due to zero point gravitons. This is implemented through an effective 
metric defined in terms of the random spacetime metric modulo a fluctuating conformal factor. 
It satisfies an effective Einstein equation coupled to an effective stress-energy tensor incorporating 
gravitational self energy of metric fluctuations as well as matter fields. The effective Einstein 
equation is expanded perturbatively up to second order nonlinearity. In the process of regularizing 
divergent integrals, a UV-cutoff is introduced whose specific value, however, does not enter into 
the resulting description of random gravity. By assuming certain physically reasonable statistical 
properties of the conformal fluctuations, it is shown that the averaged effective metric satisfies the 
empty space Einstein equation with an effective cosmological constant. This effective cosmological 
constant vanishes when only the massless matter fields are included. More generally, a finite effective 
cosmological constant compatible with the observational constraints can be obtained as long as the 
bare masses of the massive matter fields are nearly zero, or the conformal invariance of matter is 
restored at some high energy scale. 

PACS numbers: 04.20.Cv, 04.60.Bc, 05.40.-a 



I. INTRODUCTION 

Zero point energy and vacuum fluctuations have fasci- 
nated physicists 0, i, i, i, i, i] and the wider scientific 
audience @, [H, Q ever since their prediction by quantum 
field theory. Although the precise quantum nature of 
vacuum is not fully understood, vacuum fluctuations of 
the electromagnetic field have long received support from 
experiments. They include the Casimir effect (Tol [Til [l2| . 
Lamb shift of atomic energy levels [T3 , 1"l4| and sponta- 
neous emission from atoms (TH HH, Il7i. The fundamen- 
tal fluctuation-dissipation theorem [1 81 ] provides a fur- 
ther basis for the vacuum fluctuations manifested as su- 
perconductive current noise in the Joscphson junctions 
tarM [Ill- 
Even so, elements of doubt still surround the physical 
reality of the quantum vacuum. Some argue that e.g. 
the Casimir effect could be explained without the vac- 
uum energy p2| . Laboratory measurements of vacuum 
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fluctuations can at best detect differences in the vacuum 
energy by modifying their boundary conditions [23[ . It is 
generally believed that if vacuum energy is real, then it 
must gravitate in accordance with general relativity, and 
that only through the resulting gravity the net vacuum 
energy could be determined. Recently, it is shown that 
the part of vacuum energy responsible for the Casimir 
effect does indeed gravitate |2J]. However, the gravita- 
tional consequence of the total vacuum energy remains 
controversial 0, IM H3, HI • 

By Lorentz invariance, after any appropriate rcgu- 
larization, the vacuum energy can only couple to the 
Einstein equation via an effective cosmological constant. 
There is a common perception that an exceedingly large 
cosmological constant must necessarily be built up from 
all ground states up to the Planck energy density scale 
[29l ]. While a sizable cosmological constant could have 
driven the cosmic inflation in the early universe, obser- 
vations indicate its present value to be just under the 
critical density value [3(J, H3 • 

Despite the early suggested link between the vacuum 
energy of clcmentaryparticles and the cosmological con- 
stant by Zel'dovich [32j |. a detailed mechanism is lacking. 
This motivates the "dark energy" models such as the 
quintessence fields and their extensions [29|, [3l| . Leaving 
aside the original mysterious cancelation of the huge vac- 
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uum energy, these models indeed offer arguably the most 
popular current approach to the "cosmological constant 
problem" . Nevertheless, efforts to account for the ob- 
served cosmic acceleration using cosmological perturba- 
tion back-reaction without resorting to dark energy have 
continued. The need for a better understanding of non- 
linear metric perturbations [33| and quantum gravity ef- 
fects [34| is clearly highlighted by the recent progress in 
this direction. 

Spacetime metric fluctuations are an integral part of 
the quantum vacuum. However, a fundamental difficulty 
to studying these effects is the absence of a consistent 
quantum theory of gravity with an appropriate classical 
limit. Recently, progress has been made u sing a stochas- 
tic semiclassical gravity approach [H, [H, [37], HI] . How- 
ever, these results are limited to linear metric fluctuations 
and are explicitly evaluated only for the vacuum states of 
conformally invariant fields. Moffat's stochastic gravity 
[39| is nonlinear and capable of including general matter 
sources. It assumes a fluctuating gravitational coupling 
strength with certain statistical properties. 

In this paper, we address the issue of metric fluctua- 
tions and their implications on cosmic acceleration with- 
out invoking additional hypothetical fields. We shall de- 
velop a new stochastic classical gravity approach without 
modifying the Einstein equation. Rather, we seek ap- 
propriate fluctuating solutions to these equations in the 
spirit of Boyer [40l |. However, we do not advocate the 
stochastic classical approach to replace the ultimate need 
for quantum gravity. The new approach should nonethe- 
less allow us to model low energy metric fluctuations due 
to the quantum vacuum and assess their effect on clas- 
sical spacetime at a larger scale. A follow-up paper [4lJ 
will explore the possibility of decoherence effects due to 
nonlinear metric fluctuations [12, S3] ■ 

Classical stochastic behaviour is evident in many ob- 
served quantum vacuum phenomena. The effects in- 
volved may be explained using classical electromagnetic 
fields satisfying Maxwell's equations, as pointed out by 
Welton and others QJ, 53, \M, HI • The fields are how- 
ever subject to fluctuating boundary conditions, as for- 
mulated in Boyer's random electrodynamics [13, E3], to 
mimic vacuum fluctuations. York championed a novel 
gravitational analogue in which black hole entropy and 
radiance are derived from quasinormal mode metric fluc- 
tuations [il| . These fluctuations are prescribed by a clas- 
sical Vaidya geometry satisfying the Einstein equation 
with amplitudes set to the quantum zero point level. An 
improved quantum treatment of the problem using path 
integral is recently provided in [49j . 

The stochastic classical approximation to quantum 
fluctuations is physically justified for low energy vacuum 
effects. However, new subtleties arise when interactions 
and nonlinearities are included. One hopes to derive 
these effects from a full quantum theory. But, when it 
comes to gravity, there isn't one ready at hand! Given the 
classical nature of the low energy quantum fluctuations 
of matter, it seems inconsistent if they do not couple to 



gravity classically. A naive stochastic classical treatment 
of general relativity in vacuum could be applied pertur- 
batively as follows. For linearized gravity, the Boyer type 
boundary condition results in a sea of zero point gravita- 
tional waves. The effective stress-energy tensor of these 
stochastic gravitational waves then enters into the higher 
order Einstein equation [H3, HH as a source. This effec- 
tive energy density is positive and formally infinite. Even 
after a cutoff at e.g. the Planck scale, there appear to be 
problems still. The resulting effective stress-energy ten- 
sor resembles that of a massless spin-2 radiation fluid. It 
cannot be interpreted as a cosmological constant contri- 
bution, thereby breaking the Lorentz invariance of the 
resulting vacuum state. Its large energy density would 
lead either to a rapid expansion of the universe as dis- 
cussed by Weinberg [29[ or to a drastic collapse of space- 
time as noticed by Pauli [1, 0] . It seems that adding the 
vacuum energy of matter could serve only to worsen the 
situation. 

What has not been taken onboard in the argument 
above is that nonlinear metric perturbations may them- 
selves fluctuate. All experiments done so far on vacuum 
fluctuations have clearly demonstrated the crucial de- 
pendence of their physical consequences on the choice 
of boundary conditions. A physically acceptable bound- 
ary condition for the spacetime metric should allow for 
its higher order fluctuations in the presence of nonlinear 
gravitational interactions. However, cosmological obser- 
vations yield constraints on classical, averaged quantities. 
Hence the boundary conditions on the nonlinear metric 
perturbations can practically be imposed only after an 
appropriate averaging procedure. The statistical proper- 
ties of the fluctuating part of the metric are determined 
by requiring that the resulting averaged classical equa- 
tion matches the Einstein equation in empty space with 
an effective cosmological constant. 

This paper is devoted to developing and presenting 
this general framework by showing that such classical 
boundary conditions are attainable if a stochastic confor- 
mal modulation of the metric appears at the microscopic 
level. It is known from the canonical analysis of general 
relativity [13, HH that an oscillatory conformal modula- 
tion of a metric induces an effective stress-energy tensor 
with a negative kinetic term. This suggests that confor- 
mal fluctuations may yield a regularization mechanism 
for the overall amount of vacuum energy. We will show 
that, provided the conformal fluctuations satisfy certain 
statistical conditions, the vacuum structure of spacetime 
can be made compatible with Lorentz invariance at a 
classical scale together with an effective cosmological con- 
stant. 

In what follows, all calculations will be done locally 
with respect to a physical inertial laboratory frame. 
When viewed at this laboratory scale, (apparently) 
empty spacetime appears to be flat and described by a 
Minkowski metric. We shall further assume that there 
exists a cutoff scale limiting the laboratory observer ac- 
cess to the microscopic structure of spacetime and its 
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effects upon physical probes. In particular, this cut- 
off scale also fixes the energy scale which is accessible 
by the laboratory observer. Throughout this paper we 
use the signature convention (—,+,+,+) and, unless 
otherwise stated, we work in dimensionless units with 
c = G = H = 1. Moreover, when a tensor, e.g. A a b, 
is a functional of other tensors, e.g. B cc i, C e , • • ■ , we will 
sometimes use the notation A a b[B, C, • ■ ■ ], so that the ex- 
pression [B, C, ■ ■ ■ ] indicates functional dependence upon 
the tensors B cd , C e , • • • . When outside the brackets, the 
symbol B := g ab B a t, will denote the trace. With these 
conventions, the Einstein equation with a cosmological 
constant A reads G a f,[g] + Ag a b = %irT a b. 



II. MICROSCOPIC RANDOM GRAVITY 
A. Random scale and stochastic classical approach 

We consider the microscopic structure of spacctime at 
a scale I := XLp, where hp f=a 10~ 35 m is the Planck 
scale and where the dimensionless parameter A > 1. It 
sets the benchmark between the full quantum gravity 
domain and a scmiclassical domain, in which spacetime 
properties still inherit traces of the underlying quantum 
gravity physics, though being expected to be treatable 
by semiclassical means. 

To study spacetime at such a small scale we extend 
Boyer's framework of random electrodynamics S3| 
and describe spacetime vacuum fluctuations by means 
of stochastic classical fields. This kind of approach is 
necessary, as a definite quantum gravity theory is still 
lacking. Boyer's work shows that having a stochastic 
classical field results by choosing appropriate fluctuat- 
ing boundary conditions for the classical field equations. 
The classical-random field can simulate quantum vacuum 
fluctuations in the sense that a variety of physical phe- 
nomena, e.g. Casimir force between two plates, can be 
accounted for within the classical-random scenario. 

Even in empty spacctime the vacuum energy of matter 
fields is still a source of gravity. Therefore we consider 
the Einstein equation for the spacetime metric jab 

G a6 [ 7 ] = 87^^,7], (1) 

where T ab is a model stress-energy tensor describing 
the overall vacuum energy contributions coming from all 
matter fields, collectively denoted by tp. The metric "f a b 
and the matter fields ip are considered to be randomly 
fluctuating at the scale /. Accordingly we will refer to 
I as to the random scale and interpret it as the typical 
scale above which quantum vacuum properties can be 
approximately described by means of classical stochastic 
fields. Below I and closer to the Planck scale a full theory 
of quantum spacetime is required. 

Since we are considering here vacuum metric fluctua- 
tions in an otherwise empty universe, these can in prac- 
tice be assigned as random perturbations about some 



background metric g^ b . In the following calculations we 
work locally and with respect to a physical inertial lab- 
oratory frame, whose typical scale Lj a b is much larger 
than the random scale, yet small enough for the back- 
ground metric to appear Minkowski in the appropriate 
coordinate system. With this choice we have gf ib = r] ab 
and 

Jab=Vab + ^ b ) +7$ + ..., (2) 

where 7^ indicates a small fluctuating term. Here and 
henceforth we follow the standard notation where a su- 
perscript (n) denotes an n-th order perturbation. In the 
subsequent expansion of field equations, matter fields ip 
will be treated as first order quantities. 

The classical equation ([1]) will be analyzed explicitly in 
sections IIIII and IIVI for the first and second order metric 
perturbations. Boyer's type fluctuating boundary con- 
ditions will be imposed on the linearized Einstein equa- 
tion, so that the resulting randomly fluctuating tensor 
7^ can be interpreted as describing graviton fluctua- 
tions at the random scale I. Its induced stress energy 
tensor is quadratic in the fluctuations and, together with 
the part of T ab which is quadratic in the matter fields 
ip, it well enter the second order equation for the higher 
order metric perturbations. These are not expected to 
describe a physical field having its own amount of vac- 
uum fluctuations and, accordingly, they will not be linked 
to Boyer's type random boundary conditions. They will 
nonetheless be fluctuating as a result of their coupling to 
matter and graviton fluctuations. 



B. Macroscopic Lorentz invariance of vacuum 

Once the vacuum properties of spacetime at the ran- 
dom scale I are described in terms of a stochastic classical 
metric, all derived physical quantities, e.g. the Einstein 
tensor, are fluctuating also. As already noted by Boycr 
[4(i| vacuum must be Lorentz invariant when viewed at 
some appropriate macroscopic scale. As a result the only 
way it can possibly contribute to the Einstein equation is 
through an effective cosmological constant term. Lorentz 
invariance also implies that vacuum statistical properties 
must be the same regardless of space position and direc- 
tion, i.e. homogeneity and isotropy hold in a statistical 
sense. 

In order to include these properties into our formal- 
ism we want to recover classical and smooth quantities 
from the fluctuating fields. This can be obtained as a 
result of an averaging process. To this end we follow the 
spacetime averaging procedure described in [5(J, HH , so 
that fluctuating tensors average to tensors. This process 
involves a spacetime averaging over regions whose typi- 
cal dimensions are large in comparison to the fluctuations 
typical wavelengths but smaller than the scale over which 
the background geometry changes significantly. Accord- 
ingly we introduce an averaging classical scale L such 
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that I <C L <C -Liab- For completeness we also define a up- 
per breakdown scale L max > Li a b as that scale at which 
the background deviations from a flat geometry start to 
be significant. 

While keeping in mind that we are here only consider- 
ing 'apparently' empty spacetime, a final comment about 
the involved physical scales is in order. The following hi- 
erarchy holds, with Lp < I « L <C Li a b iS L max . The 
precise characterization of the classical scale L is that of 
the smallest scale at which classical, Lorentz invariant 
spacetime starts to emerge as a result of the averaging 
process. Though much larger than I, the classical scale 
L is still expected to be very small in comparison to the 
laboratory scale (table Hj). The limits of the presently 
suggested theory are thus clearly set: (i) below the ran- 
dom scale I the random fields approximation breaks down 
and a full quantum gravity theory would be needed; (ii) 
spacetime starts to be smooth and classical when viewed 
at the classical scale L; moreover, from L and up to the 
upper breakdown scale L max , the averaged fluctuations 
look like small, classical metric perturbation upon a flat 
background; (iii) beyond the scale L max , the approxima- 
tion of flat Minkowski background breaks down. 
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gravity 
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gravity 


gravity 
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TABLE I: A guide to relevant physical scales. The current formu- 
lation of the random gravity approach is supposed to yield a valid 
description of spacetime from the random scale I, up to the break- 
down scale L max . Below I a full quantum gravity theory would 
be needed. Above L max the Hat background approximation breaks 
down and one should consider the more general case of small fluc- 
tuations upon a curved geometry gfl. 



C. Conformal fluctuations and the effective 
Einstein equation 

The linearized Einstein equation with fluctuating 
boundary conditions will be connected to the vacuum 
fluctuations of graviton and will pose no problems. 
The second order equation will contain a source term, 
quadratic in the fluctuating fields ip an( l 7„j, , whose phys- 
ical effect upon the properties of empty spacetime must 
be carefully assessed. To this end, and without loss of 
generality, we introduce an effective metric g ab , confor- 
mally related to the metric "f ab via 

lab = e 2a g a b- (3) 
Hereafter we refer to a as to the conformal field. 



The equation satisfied by the effective metric g a b is 
found in a standard way from G a b[7] = 87rT a h [ip, 7] by re- 
expressing the Ricci tensor i? Q b[7] and Ricci scalar R[y] 
in terms of g ab , its compatible covariant derivative V a 
and the conformal field a as 

Rab[j] = Rab[g] - 2V a a 6 - g a bOa + 2a a a 6 - 2g ab a' c a, c 
and 

R[j] = e ~ 2a {R[g} - 6Da - 6a' c a, c } . 

Here a^ a '■= V a a = d a a, □ := V C V C , while R a b[g] and 
R[g] are the Ricci tensor and Ricci scalar for the effec- 
tive metric g ab - The Einstein tensor G a b[ 7 ] = Rab[j] — 
\R[l]"fab follows as 

Gab[l] = G ab [g]-2V a a, b + 2g ab na + 2a, a a. b + g ab a' c a, c , 

(4) 

and we have the following effective Einstein equation for 
the effective metric g ab 

G ab [g] = 8Tr{T ab + E o6 ). (5) 

The effective stress energy tensor for the conformal field 
is 

E ab := T} ab + X^ b , (6) 
where we have conveniently split H ab into the two parts: 

S afc : = J" (V a a b - gabO a) , (7) 

and 

^Ib (a. a a :b + i.9afccv' c a c ^ . (8) 

The effective stress energy tensor defined above de- 
pends on a only through its derivatives. Any physical 
situation with a = const, would yield Y> ab = 0. In this 
case the two metrics i ab and g ab are equivalent and sim- 
ply describe the same spacetime with a different choice 
of physical units. However, if a varies throughout space- 
time, the induced non vanishing stress energy tensor will 
affect the metric tensor g ab . In particular, conformal 
spacetime modulations induced by a randomly fluctuat- 
ing a would contribute, together with ip and 7^ , to the 
vacuum structure of spacetime. 

Since the metric "f ab satisfies its own Einstein equation 
C a h[7] = 8irT ab [i/j, 7], there is a fundamental arbitrari- 
ness connected to the introduction of the conformal split 
([3]). This arbitrariness has two main implications: (i) a 
complete freedom in the choice of a possible dynamical 
equation for a, whose Boyer's type boundary conditions 
would lead to a randomly fluctuating conformal field; (ii) 
once an equation is chosen the statistical amplitude of 
the conformal fluctuations would still be completely ar- 
bitrary. 
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The arbitrariness in the introduction of a will be lifted 
by imposing that the averaged second order Einstein 
equation for the metric g a b preserves vacuum Lorentz 
invariance at the classical scale and reproduces the ap- 
propriate empty space classical Einstein equation with 
an effective cosmological constant. We will show indeed 
that these requirements suggest a simple wave equation 
holding for a and fix precisely the statistical amount of 
the conformal fluctuations amplitude. 

D. Characterization of the spacetime fluctuations 

Within our framework, spacetime at the random scale 
presents two types of fluctuations: (i) conformal, scale 
fluctuations induced by the randomly fluctuating confor- 
mal field a; (ii) effective metric fluctuations as described 
by g a b- We now proceed to set up the relevant stochastic 
properties of a and g a b- 

(i) Conformal fluctuations. The field a describes small 
fluctuations in the local scale of spacetime. We assign it 
as a first order, stationary and stochastic field satisfying 

(a) = 0, (a a ) = 0, (9) 

where (•) indicates averaging. The particular choice 
(a) = corresponds to the fact that the expectation 
value of the quantum field operator a (object of the 
underlying quantum theory) in the appropriate vacuum 
state |0) should vanish. 

In section IIIII we will interpret the first order metric 
fluctuations in 7^ as describing zero point graviton fluc- 
tuations. Accordingly, we introduce the convenient no- 
tation 

@ab '■= lab ■ (10) 

It is our choice to prescribe the conformal fluctuations 
in a to be in addition to and strongly uncorrelatcd with 
the graviton fluctuations. This choice serves to simplify 
the technical derivations presented in this paper. A more 
detailed analysis shows that our main results arc never- 
theless general enough without assuming the statistical 
independence of the conformal and graviton fluctuations. 
Further comments upon this particular choice are given 
near the end of section IIIII Accordingly, we impose the 
statistical property: 

{a,..fl ab ,..) = 0, (11) 

where ", • • • " is a shorthand to indicate the derivatives. 
In particular this implies that 

(aP a b) = 0. (12) 

When /3 a b is split into its independent traceless fl* b and 
trace fl := n ab fl a b parts as 

Pab = Kb + \vabfl (13) 



the above relation also implies 

(afl* ab ) = (a/3) = 0. (14) 

Note that indeces are here raised by r\ a b as we are working 
with first order quantities. 

(ii) Effective metric fluctuations. Since empty space- 
time with vacuum fluctuations appears approximativcly 
flat at the laboratory scale Li a b, we assume a microscopic 
random effective metric of the form 

gab = Vab + g^b + 9 ab ■ ■ ■ , (15) 

where r\ a b is the Minkowski metric and g a n b a stochastic 
tensor describing n-th order fluctuations. These are small 
in the sense that, in some appropriate coordinate system, 
they satisfy \g ab \ ^ 1- We shall also use the shorthand 

qab'.= 9$>. (16) 

For each order n of the expansion, we have the classical 
(non-fluctuating) part of g a ^ given by 

9a r - ( 9a n b } ) (17) 

and define the fluctuating part of g ab by 

J«)F ._ M / (n)\ riR s 
9ab — 9ab ~ \9ab / U 8 ) 

so that g$ = g a n b )C + g a n b W . It follows that g a f F has a 
zero mean: 

(si: )F ) = o. (19) 

The average of the random effective metric at the scale 
L, yields the classical effective metric 

(9ab) =■■ 5ah = Vab + Q^b + 9aT + 0(3), (20) 

where the terms g a ^ C could account for some large scale 
deviation from flat spacetime due to vacuum fluctuations. 
In terms of the metric ^ a b we have 

lab = e 2a g a b 

= (1 + 2a + 2a 2 + . . .) x fr) ab + q a b + 9^ + ■■■)■ 
Keeping terms up to second order we have 

lab = Vab + Pab + lab + °( 3 )> ( 21 ) 

with linear fluctuations, 

flab = Qab + 2an ab (22) 
and with nonlinear metric fluctuations given by 

lab = 9ab + 2a 2 V ab + 2aq a b- (23) 
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After averaging, the corresponding classical metric 
lab :== {lab) i s > m the second order, 

lab = Vab + q°b + g a l )C + 2 (a 2 ) rj ab + 2 {aq ab ) + 0(3). 

(24) 

The mean of the cross term involving a and q a b gives 

(aq ab ) = (af3 ab ) - 2, ]ab (a 2 ) = -2r ]ab (a 2 ) , (25) 
where we have used (1121). Then 



7a C fc = (1 - 2 (a 2 )) Vab + qC b + g% C + 0(3). (26) 

Comparing with (|20|) we sec that, apart for the re-scaling 
factor 1 — 2 (a 2 ) and up to second order, the classical 
properties of the spacetime metric can be described by 
the classical perturbations q^ b and g ab of the effective 
metric g ab . The scaling factor should be positive in order 
for the metric not to become singular. Within the cur- 
rent second order approximation, this implies the condi- 
tion (a 2 ) < 1/2. This bound is however likely to be an 
artifact due to the expansion truncation, since the full 
conformal factor cxp(2a) is always positive. 

It is worth further investigating the correlation prop- 
erties between the conformal fluctuations a and the ef- 
fective metric linear perturbation q ab . By collecting all 
the traceless and trace parts in (12"2")) we have 



P*ab + \vabP = q* ab + ^r) ab (q + 8a), (27) 

where the traces are obtained as (3 = r] ab (3 ab and q = 
ij ab( lab- The above equation implies the two algebraic 
constraints 



€b = Kb 

q = /3-8a. 



(28) 



The first constraint involving the traceless parts implies 



(29) 



as we know from (fT4|) that a and (3* b are uncorrelated. 
From this property and equation (|25|) we get 



(aq) 



.(a 2 ) 



(30) 



The last two equations are important in that they show 
that the conformal field is really just correlated to the 
trace part of the linear metric perturbation q ab . 

We conclude this section by remarking that, if the met- 
ric perturbation (3 a b was put into a TT gauge, then (3 = 
identically and, as a consequence, q = —8a. In this case, 
apart for a multiplicative factor, a would coincide with 
the trace of q ab . It follows that, once (3 = 0, q ab cannot 
be put into a TT gauge as well, or it would be a = 0, 
leading to the trivial case g ab = "f ab . Considering that 
gravitons would be described in the TT gauge by (3* ab and 



also in view of the identity q* b = (3* b , we can interpret the 
traceless part of the metric fluctuation q ab as describing 
graviton in the spacetime (Ai, g a b)- We also remark that 
the assignment of a as an external fluctuating field effec- 
tively attaches physical meaning to the trace fluctuations 
of q a b- 



E. Matter field stress-energy tensors 

Vacuum carries contributions from all physical fields of 
nature. As a result the random gravity framework is not 
complete without inclusion of the matter fields. These 
are considered to be in their ground state. We remark 
that the vacuum state for the various matter fields is well 
defined on a flat Minkowski background. However, since 
we work at the random scale I, we represent matter fields 
vacuum fluctuations by means of first order and stochas- 
tic quantities. At this stage we adopt a minimum ap- 
proach by neglecting the non-gravitational interactions 
between various matter fields components. The stress 
energy tensor T ab in ^ describes matter fields and, ide- 
ally, carries contributions from all sectors of the Standard 
Model. Then T ab = • T^ b , where the index j runs over 
all matter fields. 

The detailed microscopic expression of the generic 
component T 3 ab will depend quadratically upon the corre- 
sponding fluctuating matter field, as well as on the ran- 
dom metric jab = Vab + J2 n lot? ■ ^ s a rGSU lt the stress 
energy tensor at the random scale I is a also a stochas- 
tic quantity. The dependence on the 7^ would account 
for the back-reaction of gravity fluctuations upon matter 
fields. However, as long as we work up to second or- 
der, only the flat classic background r] ab will appear but 
not the metric fluctuations 7^. In this sense it is like 
having fluctuating fields on a Minkowski background and 
the effect of gravity upon the stress energy tensor would 
only appear as a third order effect. Within this frame- 
work, in an appropriate coordinate system, matter fields 
stress energy tensor component are quadratic fluctuating 
quantities defined on a flat background. 

The corresponding energy density contribution can be 
defined at the classical scale L in a statistical sense 
through an averaging procedure. Provided the high 
frequency components are cutoff at the random scale 
I = XLp, the average will be well defined and finite. 
Then the quantity (T ab ) gives rise to a macroscopic stress- 
energy tensor at the classical scale L. Because of homo- 
geneity and isotropy holding at the classical scale, its 
most general form is 



(T ab 



( p 
Op 







p 

V p J 



(31) 



where p := (Too) is the energy density and p := | (T%), 
i = 1,2,3 is the pressure. The dominant energy condi- 
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tion, i.e. p > and p > p, is normally thought to be valid 
for all known reasonable forms of matter [55| , at least as 
long as the adiabatic speed of sound dp /dp is less than 
the speed of light. This is true for massless fields since, 
in this case, p = 3p. In appendix [A] we show that the 
stronger condition p > 3p is satisfied by massive fields in 
their vacuum state, at least in the ideal case in which in- 
teractions can be neglected and the field masses are much 
smaller than the Planck mass. 



F. From random gravity to classical gravity 

The physical picture is such that, at the random scale 
/ = XLp, spacetimc and matter fields arc randomly fluc- 
tuating. As a result we have the random Einstein equa- 
tion for the effective metric g ab 

G ab \g]=^T^,a,g], (32) 

with ip denoting collectively all matter fields and where 
the effective stress energy tensor induced by the confor- 
mal fluctuations is 

TabbP, a :9] = SafeKs] + Z 2 ab [a,g] +T ab [ip,a,g], (33) 

the linear and quadratic parts in a being given in equa- 
tions (Zj) and ©. 

At the classical scale L and above, the averaged mi- 
croscopic equation will yield a corresponding classical 
equation in terms of smooth, non-random quantities 
(G a b[g]) = 87r (T^[ip, a, g]) . We will now proceed to ap- 
ply an order by order expansion scheme to equation (|32p 
and study under which conditions it takes the only form 
which is compatible with Lorentz invariance of vacuum, 
i.e. that of the empty space classical Einstein equation 
with an effective cosmological constant term. 

III. FIRST ORDER ANALYSIS OF THE 
RANDOM EINSTEIN EQUATION 

A. Linear solution and graviton vacuum 
fluctuations 

We now proceed to linearize the microscopic random 
Einstein equation (|3"2")) and study the structure of its ef- 
fective stress energy tensor, including vacuum conformal 
fluctuations. The terms T? ab and T ab contain fluctuations 
starting from second order so that only the term E^ b 
contributes. Using equation J7J) we have the linearized 
random Einstein equation 

G { ab [?] = 2d a9 b a - 2 Vab d c d c a, (34) 

where G~) is the usual linear operator resulting from the 
linearized Einstein tensor (5ll | and q ab represents the first 
order metric fluctuations. Here the indices are shifted 
using rj ab and r/ ab . 



This equation can be simplified. Indeed, from q ab = 
f3 ab — 2ar] abl it is readily verified using the explicit ex- 
pression for G^b given in appendix IB1 that 

G { ab M = G { a b ] \P - 2a V ] = G§ [(3} + 2d a d b a - 2 Vab d c d c a. 

(35) 

Comparing with (|34[) we see that the first order approx- 
imation to the microscopic Einstein equation takes the 
form 

Gi 1 5 ) [/?] = 0, (36) 
with (3 ab = q ab + 2ar/ ab and 

{(3 ab ) = (q ab ) = q%. (37) 

In the Lorentz gauge, equation (|3G[) is the ordinary 
wave equation and it is usually considered to describe 
weak gravitational waves (GWs). Within our current 
random framework, it can be assigned fluctuating bound- 
ary conditions. As a result, the fluctuating tensor /3 ab is 
thought to represent the vacuum fluctuations of gravi- 
ton. Note that, as explained in section Hi C[ the confor- 
mal field a can in principle be assigned arbitrarily. Once 
this is done, q ab is fixed from q ab = /3 ab — 2arj ab in such 
a way that |34|) is automatically satisfied. We remark 
how the conformal fluctuations act in equation (|34|) as 
a fluctuating forcing source term, implying that the first 
order effective metric fluctuations q ab correlate to a. This 
observation rules out the possibility (aq ab ) = and mo- 
tivates our choice (a/3 ab } = in section TlID I as. in fact, a 
physically reasonable scenario. 

The corresponding equation holding at the classical 
scale L is found by taking the mean in ([34|) or |36|) . Since 
G^a b is a linear operator it commutes with the average 
operation [U and {G§\0\) = G$M], implying 

G^[q C ] = 0, (38) 

where we have used (|37|) . The first order, classical, met- 
ric correction q^ b would describe classical gravitational 
waves propagating on a Minkowski background. Since we 
want to characterize spacetimc vacuum properties only, 
we set these GWs to zero by choosing 

q^b = o. (39) 

This represents the first order solution to our problem. 
We have here the important result that, in the first order, 
vacuum fluctuations do not give any visible effect on the 
structure of spacetimc. 

B. Drift of classical spacetime due to non linear 
vacuum effects 

After taking into account the first order solution, the 
metric structure follows from equations (|2l"j) - (f2"3"]) as 

lab = r]ab + ^6 + 7^+0(3), (40) 
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with the zero mean, linear fluctuations, 

Plb = Qab + (41) 

describing graviton vacuum fluctuations, and with second 
order fluctuations 

lab = 9 { ab + 2 " 2 'M + 2ag a F fc . (42) 
After averaging, the full classical metric follows as 

< 7ab ) = (1 - 2 (a 2 )) Vab + g$ C + 0(3), (43) 

where we have used (|25[) . This result shows that a de- 
viation from a flat Minkowski background can arise as a 
second order drift effect due to nonlinear vacuum fluctu- 
ations. The classical equation satisfied by g^J C is found 
by a second order analysis of the averaged random Ein- 
stein equation. 

IV. SECOND ORDER ANALYSIS OF THE 
RANDOM EINSTEIN EQUATION 

A. Second order equation and gravitons effective 
stress energy tensor in vacuum 

Expanding the fluctuating Einstein tensor in the ran- 
dom equation (|3"2")l up to second order yields the result 

m 

G ab = G$ [q] + G$ [g^] + G<£ [q] + 0(3), (44) 

where G^ b is the linear operator introduced above while 

G^ab ^ s quadratic in its argument [56[ and its explicit 
structure is reported in appendix [Bj 

Because observations of the spacetime metric involve 
classical, macroscopic scales, we want to analyze the 
structure of the averaged second order equation. Taking 
the mean in Q44|) , the first term vanishes by virtue of the 
first order solution. Then, using the fact that q ab = q^ b 
we get 

(Gab) = G§\g<» c ] + (G2V]) + 0(3). (45) 
Using ql b = /3 F b - 2arj ab we obtain 
(Gab) - Gl%^} + (G^)) + (G»N) . (46) 

A fourth term deriving from (GQ [q F ]) and containing 
products of the derivatives of /3 F b and a averages to zero 
because of the fact that these two quantities are strongly 
uncorrelated as expressed in 

The second term in (|46p gives rise to the effective 
stress-energy tensor 

(T%»[f3 F ]):=-±{G$^}), (47) 



which is quadratic in the first order fluctuations /3 F b and 
acts, together with matter and conformal field fluctua- 
tions effective tensors, as an extra source for the classical, 
second order metric drift term g^ b ■ 

The structure of this effective stress-energy tensor has 
been studied extensively in [5(j, where it is shown to be 
positive definite, traceless, and well defined in describing 
the energy content of linear gravitational fluctuations, as 
long as the wavelengths of the fluctuations involved are 
shorter than the typical scale over which the geometry 
of the background metric varies significantly. In our case 
this means that the graviton fluctuations wavelength con- 
tained in /3 F b should be shorter than the breakdown scale 
L mffl . This would technically imply an infrared cutoff. 
However, since in the current version of the theory we 
are working with a flat background we will ignore this. 
Within our present framework, (T^ b w ') will be connected 
to the energy density due to graviton vacuum fluctua- 
tions. Given the properties of (T^ b w ^ and the fact that 
it effectively describes gravitons as massless spin-2 parti- 
cles, we can avoid considering its explicit averaged struc- 
ture by simply including gravitons within the collection 
of matter fields, as described collectively by tjj. 

With this in mind, the averaged second order equation 
takes the form 

GilV 2)C ] = 8n(T^[a,f] (48) 

where tp now includes gravitons as represented by /3 F b 
and the superscripts ^ indicate second order quantities. 
The first term on the r.h.s. is the second order, aver- 
aged effective stress-energy tensor due to the conformal 
fluctuations, given by 

- (h G< $ [ar,] ) ■ (49) 

B. Effective stress energy tensor analysis 

1. Conformal field 

We analyze the three contributions to (T^b ) given 
in (|4"5|) . More details on the derivations of the following 
equations are given in appendix [Bj From the first 
term can be seen to contain second order terms involving 
products of a and g F b . Using equation (|4"Tj) to replace q F b 
and using the fact that a and /3 F b are uncorrelated, this 
term turns out to be independent of /3 F b and takes the 
form: 

(x£< a) [a, /3 F ]) = i- (2a a a b + 77 afc a< c a c } . (50) 

The second term in (|49l) follows from (jSJ) replacing g ab 
with riab' 

( S afc 2) H) = (<x, a a, b + ^r]aba' c a c j ■ (51) 
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(2) 

Using the explicit form of the non linear operator G ab , 
the third term can be calculated to be: 



1 



2tt \ a0 L '7 4tt 
Collecting the above results, equation (|4"9")l becomes 



(52) 



T CF(2) 



i 



l 



aj ; = — — \ a, a a,6 - -VabW a,, 



(53) 



Considering the stochastic properties of the confor- 
mal fluctuations in different spacetime directions to be 
uncorrelated and since (a_ c ) = we have (a i0 Q j) = 
(a, a ) (o.b) = <^ab ( a2 a) i where <5 a b is the four dimensional 
Kroneckcr symbol. Isotropy at the classical scale implies 
that, for all three space directions, 



(a 2 ) = (|Va| 2 )/3, i = 1,2,3, 



(54) 



where (| Va| 2 ) = (a 2 ) + (a 2 .) + (a|) . Then it is readily 
verified that the average of the tensor in (|53p has the 
structure 



where 



I Pa 




V o 



<« o) + 






3 / 



(55) 



<|Va| 2 )' 



< 



(56) 



Quite remarkably it is traceless 



0. 



The classical, average quantity p a corresponds to an ef- 
fective, negative energy density connected to the vacuum 
conformal fluctuations. 



2. Matter fields 

The second order, averaged matter stress-energy tensor 
is given in equation (|3ip in terms of the matter vacuum 
energy density p and pressure p. As already discussed, 
it now also includes gravitons as massless particles. Per- 
forming a trace decomposition it has the structure 






V o 





3 







3 





\ 




p± J 

3 ' 



+ \v ab (T^) (58) 



where 



P% : = ^ (P + P) ■ 



(59) 



is the energy density of the traceless part and 

T^='3p-p. (60) 

is the trace part, with = 7? a& T a (2) . 

To proceed, it is conceptually useful to separate the 
fields ijj into massless fields (including gravitons) and 
massive fields, collectively denoted by '5 and T respec- 
tively. Massless fields contained in if? must have a trace- 
less stress-energy tensor. This implies the usual equation 
of state 



p<s, = 3p*, 



(61) 



where p* and py denote the contribution to the vac- 
uum energy density and pressure coming form all mass- 
less fields, including graviton. In a similar way px and 
px will indicate the corresponding contributions from the 
massive fields. With this notation equations ([59| and 
dnOJ) can be re-expressed as 



p% := p* + - (pr + Pr) > 0, 



T (2 



3pi 



Pr, 



(62) 



(63) 



showing explicitly that only the massive fields contribute 
to the trace part of matter vacuum energy density. 
Within our level of approximation, trace anomalies due 
to either matter interactions or metric fluctuations are 
not taken into account. 

Finally, the averaged equation holding at the classi- 
cal scale L and fixing the classical second order metric 



( 2)C 

(57) perturbation g ab is 



G^[ 9 ^ G }=S^(T^[a]+T^,T] 



(64) 



where (T^ b F ^ [a]) is given in 
" and ([S3]). 



and (7i 2) [*,T]) 



C. Vacuum energy balance equation 

The classical, second order equation (f6"4")> with no con- 
formal fluctuations (a = const. =>• S a b = 0) would lead 



to G 



,^[gW c ) = 8WT^[tf,T]V Given the structure 

of the vacuum, averaged matter stress energy tensor in 
(|58[) we immediately notice a problem connected to its 
traceless part. Indeed, even with the adopted ultraviolet 
cutoff at the random scale I, the energy density p*^ will 
attain a huge value and, more interestingly, will break the 
Lorentz invariance of vacuum at the classical scale. Even 
in the idealized situation where the matter fields were to 
be neglected, the traceless effective stress energy tensor 
associated with graviton, included in would still cause 
the same problem. 
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The above conclusion can be avoided once the effects 
of the nonlinear metric fluctuations are carefully taken 
into account. When these are prescribed to include some 
amount of conformal fluctuations, Lorentz invariance can 
be restored if the following vacuum energy balance equa- 
tion holds 



= 0. 



(65) 



Then the only residual vacuum effect at the classical scale 
would come in the form of an effective cosmological con- 
stant Am as defined by 



A 



M 



_ 27r ^(r( 2 )[T] 



8ttG 



-PM, 



(66) 



where the matter related quantity pm has the dimensions 
of an energy density and is defined as 



Pm ■= ^(px - 3p T ) 



(67) 



For clarity we have explicitly displayed in (f6"6"|) the factors 
that convert a cosmological constant into the correspond- 
ing energy density. 

We will study the implications of equations (f65|) and 
(|6"6"|) in the following sections. They embody the choice of 
the classical boundary conditions for the classical metric 
perturbation gr fe m the sense that, once they can be 
shown to hold, then the classical second order Einstein 
equation (|64|) in empty space reads 



(68) 



This formally corresponds to the linear approximation of 
classical Einstein equation with a cosmological constant 
A and an approximatively flat metric. Indeed, if g a b ~ 
Vab + h a b, with \h a b\ <C 1, then the linear terms of G a b = 
—Ag a b would yield G^[h] = —Ar) a b- This has precisely 
the same form as in (|68p , provided that we identify h a b = 



(2)C 



and A = A 



M ■ 



then (|69|) implies the following observational constraint 
on the present epoch cosmological constant 



< A obs < 10~ 125 Aj 



(71) 



The upper bound is extremely small, especially in com- 
parison to the Quantum Field Theory (QFT) theoretical 
expectations usually reported in the literature. There 
it is commonly argued that, applying QFT within the 
Standard Model and up to the Planck scale, the expected 
cosmological constant would result to be of the order of 
Ap, hence the famous discrepancy of about 120 orders of 
magnitude with the observed value (3l| . 

Our random gravity framework also predicts that an 
effective cosmological constant should indeed emerge due 
to vacuum energy effects. However, at odds with the 
usual claims, only the massive fields are expected to give 
a contribution. A necessary requirement for the observa- 
tional constraint not to be violated is that Am be posi- 
tive. This implies that massive matter fields in the vac- 
uum state should satisfy the following condition 

PM = ^- = \(pr-3p r )>0. (72) 

In appendix [X] it is shown that, in a first approxima- 
tion that treats matter fields in their vacuum state as 
a collection of non-interacting harmonic oscillators, pyi 
is positive and proportional to a squared mass parame- 
ter M 2 , defined in equation (|A22[) as a weighted average 
of the squared masses of all particles comprised in T. 
Specifically, we obtain in (|A24j) and (|A25|) that 



Pm 



A m = 



N T p P ( M 
8A 2 VMp" 



N T A P ( M 
\~Mp 



8A 2 



(73) 



(74) 



D. Effective cosmological constant 

The current estimated amount of observed vacuum en- 
ergy in the present epoch universe is [3l| 



obs 
r vacuum 



5 x 10" 



10 



-125 



PP, 



(69) 



where pp m 10 19 GeV /Lip is the Planck energy density. 
This value is deduced from the observed large scale prop- 
erties of the universe, including the approximatively flat 
geometry, as suggested by CMB anisotropics [57| . and 
the present, small, cosmic acceleration indicated by type 
la supernovae observations [EH ■ Defining the Planck cos- 
mological constant as 



A, 



8ttG 



(70) 



The above clearly shows that pm oc - 
M P « 10 19 GeV is the Planck mass 



(w?) . where 
This value can 
still be large and formally diverges as A decreases. If 
a Planck cutoff scale A ~ 1 were to be applied applied 
as per the "usual practice" , we would end up with an 
effective cosmological constant Am much smaller than 
the Planck value Ap, roughly by a factor of (M/Mp) 2 . 
Estimating this using the maximum mass scale of the 
Standard Model, namely M = M H w 100 GeV, then the 
reduction factor is (Afy/Mp) 2 ~ 10~ 34 and the resulting 
actual value for Am still would exceed the observational 
value A obs by about 90 orders of magnitude. However our 
present result, indicating that the cosmological constant 
seems to be related to the matter fields masses, opens up 
new a prospect for considerations that will be discussed 
in the final section of this paper. 
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V. ANALYSIS OF THE BALANCE EQUATION 

A. Lorentz invariance of vacuum and conformal 
fluctuations 

Turning now to the issue of Lorentz invariance, we can 
use (|56|) and the balance equation (|65|) to get 

3(a 2 ) + (|V a | 2 )=87r P ;. (75) 

In principle, this equation can always be satisfied and 
comes in the form of a statistical condition for the con- 
formal fluctuations a. It requires that their "intensity" 
be precisely adjusted so to balance the traceless part of 
the vacuum energy contribution due to all the other fluc- 
tuating fields. At the present stage this requirement is 
imposed "from the outside" , on the basis that Lorentz 
invariance of vacuum should hold at the classical scale. 
However we would expect that a high energy, underlying 
quantum theory of gravity that incorporates conformal 
fluctuations as an essential ingredient, could lead to a 
more detailed and fundamental "internal" explanation 
for the balancing mechanism to hold. Some preliminary 
but promising indications along these lines are given in 
[531 I59I I . where the canonical formulation of gravity is 
enriched to include conformal symmetry in both the ge- 
ometrodynamics and connection variables approaches. 

B. Characterizing a using the wave equation 

The conformal metric fluctuations provide a mecha- 
nism to restore Lorentz invariance of vacuum at the clas- 
sical scale provided they have suitable statistical proper- 
ties as expressed by the averages (a 2 ) and (|Va| 2 ). We 
recall at this point that since the conformal field a ap- 
pears in our present formalism as an external field that 
we can prescribe at will, we are free to assign it in such a 
way that it satisfies an extra dynamical constraint with- 
out affecting the Einstein equation. To be guided in our 
choice we notice that the second order effective stress en- 
ergy tensor in (|53"|) is traceless and depends solely on the 
first derivatives of a. If it had to have a trace term, the 
simplest and most natural combination containing first 
derivatives of a that would serve the purpose is 

a' c a c = -a 2 + |Va| 2 . (76) 

Therefore a simple and reasonable statistical condition 
to impose on the conformal fluctuations appears quite 
naturally to be 

(a- c a c ) = - (a 2 ) + (|Va| 2 ) = 0. (77) 

Then a natural choice for the extra dynamical con- 
straint is suggested by the fact that (i) a is a scalar field 
and (ii) it is expected to satisfy the statistical relation 
C3- It is 

d c d c a = 0, (78) 



i.e. the simple wave equation on Minkowski spacetimc for 
a massless scalar field. Indeed it is readily verified that, 
provided the statistical properties of a are stationary, 
equation (|75|) guarantees the condition (a — (|Va| 2 ) = 
to hold. Our choice is further motivated by the fact that 
the usual stress-energy tensor for a massless scalar field 
(f> on a flat spacetime has the form 

TT ■= ~ J W^cJ + \vab (~^c) • (79) 

Provided (4> c 4> c ) = and modulo the irrelevant rescaling 
factor of 1/27T, its average matches precisely the expres- 
sion in (|53|) . but with the negative sign. 

Once the wave equation for a is introduced, it can be 
solved with Boyer's type (random) boundary condition 
and the overall (classical) statistical amplitudes set in 
such a way that the balance equation is satisfied. It is 
precisely in this sense that the conformal random fluctu- 
ations are independent of graviton's and those of other 
matter fields. However, their overall amount is fixed at 
the classical level by the requirement of Lorentz invari- 
ance alone. Using ([77| the balance equation takes the 
very simple form 

(|Va| 2 ) = 27r P ;. (80) 

An important physical justification of the wave equa- 
tion (|78p for alpha is that its solution admits a Lorentz 
invariant spectral density suitable for describing vacuum 
fluctuations. Further, the form of this spectral density 
allows the effective stress-energy tensor of the conformal 
field to regularize the vacuum energy of matter as de- 
tailed below. 

The formal analogy between the stress-energy tensors 
([79")) and ([53")) and equation ([50)) above suggest that we 
introduce a rescaled conformal field as 

V '=W* (81) 

With this notation the balance equation reads 

(|V^| 2 )=p;. (82) 

C. Isotropic power spectrum of the conformal 
fluctuations 

The wave equation implies the usual dispersion rela- 
tion lu 2 = k 2 relating frequency u> and wave number 
magnitude k := |k|. The conformal field can accord- 
ingly be expressed in Fourier components whose corre- 
sponding power spectral density S v (k) can be defined in 
a statistical sense {4lj . Since vacuum is isotropic at the 
classical scale we have accordingly S v = S v (k). It is then 
straightforward to show that 

(9 2 ) = -Jrs Jd 3 kS v (k), (83) 
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<|V^| 2 ) 



(2tt) s 



(2tt) 3 



d 3 kk 2 S v (k), (84) 



d 3 kuj 2 (k)S ip (k), (85) 



which, given the dispersion relation, are of course com- 
patible with ([77)1 . 

Since j d 3 k = 4ir f dkk 2 these imply an integration 
over the wavelengths of the conformal fluctuations. In the 
present context we shall adopt a UV, i.e. high frequency, 
cutoff. We take the corresponding minimum wavelength 
to be equal to the random scale / = XL p. This is because 
our random formalism based on the concept of classical 
fluctuating fields breaks down for scales shorter that I. 
We thus have the wave number cutoff 



kx := 



2tt 
XL~p 



(86) 



From the physical view point the dimensionless cutoff pa- 
rameter A should be in the range of A > 10°, below which 
the quantum spacetimc effects arc expected to be impor- 
tant. However, as far as Lorentz invariance of vacuum 
is concerned, it can also be thought as a regularization 
parameter whose precise value is not critical. This is so 
because the same UV cutoff is to be applied to the cal- 
culation of the matter related energy density and pres- 
sure. As a result, although the resulting expressions for 
(|V<y9| 2 ) and formally diverge as A — > 0, their ratio 
and the related balancing mechanism turn out to be cut- 
off independent. 

The statistical quantities defined above must be invari- 
ant under Lorentz transformations. To this end we take 



S v (k) :-- 



S Q hG 1 



(87) 



where, being ip dimensionless, the combination of h, G 
and c gives the right dimensions for a power spectrum 
(i.e. L 3 ), and So is a dimensionless constant that con- 
trols the overall 'normalization'. Our choice is motivated 
by the fact that this yields the Lorentz invariant mea- 
sure d 3 k/u)(k) [60l ] and it implies an energy spectrum 
p(lj) tx w 3 , which is shown by Boyer (see ref. below) to 
be the only Lorentz invariant vacuum energy spectrum 
for a massless field. Of course Lorentz invariance is pre- 
served provided the cutoff k\ is given by the same number 
for all inertial observers, as discussed in details by Boyer 



effective transition scale from the random to the purely 
quantum domain. In principle, this effective cutoff value 
could be estimated, for example, through high preci- 
sion measurement of the decoherence suffered by massive 
quantum particles propagating in vacuum [i^, IH, [lil, H^] . 
This possibility will be further investigated in the follow 
up paper 



41] 



Restoring full physical units by inserting back the con- 
stants c, h and G, equation ([82]) is 



Pi, 



p P L P {\\7<p\ 2 ) 



where Lp is the Planck length. The dispersion relation 
with restored physical constants is uj(k) = ck and the 
integrals ([Ml)- (El) yield 



W 2 ) 



So 

A 2 ' 



/, |2 \ 27T 2 So 

Substituting ([50)1 into ([88]) we have immediately 



(89) 



(90) 



(91) 



As discussed above, the cutoff parameter A dependence 
is only apparent. The matter fields "traceless" energy 
density pi is calculated in appendix [Al as 



n 2 Np P 



(92) 



where N is the total number of independent matter fields 
modes. Substituting into ([9"T[) . this leads to the cutoff 
independent result 



So = 



N 



(93) 



representing precisely the power spectrum normalization 
that is needed to balance matter fields, traceless vacuum 
energy density as described by p*^ and restore Lorentz 
invariance. 

Plugging this value back into equation (f89| . the av- 
erage conformal fluctuations squared amplitude follows 
as 



N 

2A 2 ' 



(94) 



D. Explicit solution of the balance equation 

Thanks to equations ([8"3"[) - ([8"5[) all the averages involv- 
ing the conformal fluctuations simply depend upon the 
spectral parameter So and, of course, the cutoff parame- 
ter A. Although it functions here as a regularization pa- 
rameter, we recall that a specific value for A may mark an 



VI. INTERPRETATION AND DISCUSSION 

If vacuum energy has any physical reality then it 
should act as a source of gravity according to the Ein- 
stein equation. Quantum field theory predicts a large 
amount of vacuum energy that may lead to a Planck size 
cosmological constant. However current observations of a 
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small cosmic acceleration suggest that the overall amount 
of vacuum energy in the universe must be many order of 
magnitudes smaller than predicted by the standard ap- 
proach. 

We have provided a general random gravity framework 
in which spacctime vacuum fluctuations are described by 
means of a stochastic metric satisfying the ordinary Ein- 
stein field equation. As it was shown by Boyer for the 
electromagnetic field, this is equivalent to specifying ap- 
propriate fluctuating boundary conditions. The result- 
ing fluctuating classical fields mimic quantum fluctua- 
tions, in the sense that a variety of related phenomena 
can be accounted for. In this work we assume this to be 
valid in the case of spacetime fluctuations, up to some 
energy scale which is set by the random scale close to the 
Planck scale. These are described by fluctuating metric 
perturbations defined on a flat background metric. The 
resulting random gravity theory is built so as to pro- 
vide a description of the apparent empty spacetime from 
the random scale up to some large scale L max , where 
the deviations of the background from a flat Minkowski 
start to be significant. To this end, a spacetime averag- 
ing procedure is employed to describe the passage from 
the fluctuating metric to a classical metric that describes 
spacetime beyond some classical scale L^> I. 

It turns out that, in order for the resulting theory to 
be compatible with the Lorentz invariance of vacuum at 
the classical scale, the fluctuating metric must execute 
conformal fluctuations spontaneously. This yields a met- 
ric of the form "f a b — cxp(2a)g a b- This equation must 
properly be interpreted as a simple statement about the 
two metric tensors j a b and g a b being conformally related. 
The issue of which of them actually "contains" the con- 
formal fluctuations is an interesting one. In the present 
framework this questions is partially answered by pre- 
scribing the fluctuations in a to be independent from 
the fluctuations in the first order metric perturbation 
Pab = lab j supposed to mimic graviton vacuum fluc- 
tuations. Once this is done, a appears to be related to 
the trace of q a b = g^ b \ i.e. the first order perturbation 
of the effective metric g a b- In principle we can cast (3 a b 
into a TT gauge, so that /3 = and a = — q/8, where 
f3 and q denote the traces. In this sense we can claim 
that conformal fluctuations are actually "contained" , to 
first order, within the effective metric g a b but not within 
the metric "f a b- However, conformal fluctuations prop- 
agate through nonlinearly coupling in such a way that 
a should affect the higher order metric perturbations of 
both metric tensors g a b and 7 a h. 

In the first order approximation a does not affect gravi- 
ton fluctuations and it just acts as a linear source of ran- 
dom noise for q a b- This happens in such a way that, at 
the classical scale, the classical tensors (f3 a b) and (q a b) co- 
incide and satisfy the usual linearized equation for GWs 
with no extra sources. In the second order the conformal 
fluctuations described by a induce a traceless, negative 
definite, effective stress energy tensor which acts as a 
source for the effective metric g a b- Under some circum- 



stances, this can counterbalance the predominant part 
of vacuum energy due to matter fields and coming from 
the traceless part of their associated vacuum stress en- 
ergy tensor. Without such a counterbalancing mecha- 
nism, the associated stress-energy tensor attains a large, 
close to Planckian magnitude, thereby breaking Lorentz 
invariance, even with an UV cutoff set at the random 
scale / = XL p. 

We showed that Lorentz invariance can be restored by 
the conformal fluctuations by means of a simple balanc- 
ing mechanism, in such a way that the large, traceless, 
matter fields contribution is canceled out by the negative, 
conformal contribution. This regularization process fixes 
the averaged conformal fluctuation squared amplitude 
unambiguously once these are chosen to have a Lorentz 
invariant power spectrum and to satisfy certain specific 
statistical properties. These are shown to be compatible 
with a simple massless wave equation ruling the confor- 
mal fluctuations dynamics. This extra dynamical con- 
straint is consistent with the Einstein equation since the 
conformal field is originally prescribed as an external field 
and, in principle, can be assigned arbitrarily. The main 
choice related to its statistical independence with gravi- 
ton fluctuations at the random level, together with the 
Lorentz invariance requirement at the classical level, ef- 
fectively attach physical meaning to the conformal field 
fluctuations and fix their statistical properties. 

The dimcnsionlcss cutoff parameter A introduced in the 
process is likely to mark the transition from the random 
domain to the purely quantum domain at the Planck 
scale. In this sense it sets the small scale breakdown 
limit to our current theoretical approach. However it is 
very interesting to note that, in a similar way as to what 
happens with the Casimir effect, the spectral density of 
the conformal fluctuations that results from the balance 
mechanism is cutoff independent. This is a strong result, 
"almost" entirely based upon general relativity alone. 
The "almost" simply referring to the fact that, to im- 
plement the regularization process, it has been necessary 
to modify the classical boundary conditions and enter 
into the random fields domain. 

The final outcome of the regularization process is (i) 
a Lorentz invariant "empty" spacetime at all classical 
scales and (ii) an effective cosmological constant Am, 
emerging as a second order nonlinear effect due to the 
trace part of the matter fields vacuum stress-energy ten- 
sor. The effective cosmological constant thus appears as 
the only remnant effect, built up from the quantum vac- 
uum, which could in principle have an effect at a large 
scale. At odds with what commonly stated, the mass- 
less fields do not seem to play any role in contributing 
to the effective cosmological constant, at least within the 
present second order nonlinear approach. Our final equa- 
tion for Am shows that it is proportional to a squared 
mass parameter M 2 , obtained as a weighted average over 
all massive fields. 

Even when a usual cutoff at the Planck scale (i.e. 
A ~ 1) is applied, our formula yields a leftover, effec- 
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tive cosmological constant which is roughly 30 orders of 
magnitude lower than the usual Planckian value, though 
still much larger than the almost vanishing value value 
A obs ps (8ttG/c 4 )10~ 125 p P suggested by cosmological ob- 
servations. Our current result could however be rein- 
terpreted under a slightly different and interesting per- 
spective: by identifying Am = A obs , we get an equa- 
tion for the matter field averaged effective mass M\ — 
M/(M P /\), namely 



^obs 



n ohs 

r vacuum 



(Ap/A*) VOVA 4 )' 



(95) 



Once the n umbers are plugged in this yields the result 
M ps AMpVlO" 125 ps A x 1(T 44 GeV. This result shows 
that, for any realistic value of the cutoff parameter, in- 
cluding the Planckian value Awl, the average mat- 
ter fields effective mass matching the currently observed 
cosmological constant value is basically vanishing. This 
result is interesting since it links to some scenarios dis- 
cussed in the literature and suggesting that, at the most 
fundamental level, matter fields bare masses may in fact 
be zero [621, 63] . The observed particle masses would ap- 
pear as a result of interactions. One example would be a 
Higgs, complex scalar field based interactio n 164 . [65l . 1661 ] . 
but other possibilities are being considered 63, 67l |. 

The particular approach used in this paper, where mat- 
ter interactions are not included, could be relevant in 
this perspective if we interpret the effective mass M as 
the average bare mass of the matter fields in their vac- 
uum state. If this is zero as suggested by some authors 
and, we add, if the actual mass generating mechanisms 
are not effective in the vacuum state, then a possible 
solution to the riddle of the cosmological constant prob- 
lem could be "almost" at hand. We say "almost" since 
any sound statement would first require a solid and deep 
understanding of the nature of vacuum and, especially, 
of the meaning of mass and how it comes to be. At 
this stage we point out however that, with the conformal 
fluctuations as a possible mechanism to yield Lorcntz in- 
variance and neutralize most part of the matter fields 
vacuum energy density, the possible connection between 
cosmological constant and effective (bare) matter fields 
mass is very interesting. 

A widely accepted present view is this: While physics 
may well restore conformal invariance at a unified high 
energy scale where all fundamental particles become 
masslcss, these particles nonetheless gain an effective 
mass by settling into a local minimum of the Higgs po- 
tential (or any other mass generating mechanism) at low 
energies. This implies that the effective cosmological con- 
stant Am is built from nonzero mass up to the Higgs (or 
equivalent) energy scale, say Mh ~ lOOGeV, correspond- 
ing to an effective cutoff at A ~ M P /M H in ([71)). This 
equation, together with M ~ Mh then yields the esti- 
mated effective cosmological constant to be 



which still exceeds the observational value by about 57 
orders of magnitude. However, this estimate assumes a 
simplistic zero value of the local minimum of the Higgs 
potential. It has been suggested that, it is the local max- 
imum, where the Higgs field vanishes, that should be 
zero [20]. While other "calibrations" may exist, it is in- 
teresting to note that if this is true, then the equilib- 
rium potential energy of the Higgs field would amount 
to an additional effective cosmological constant, that is 
(roughly within the same order of magnitude) equal and 
opposite to Am given in ([M]) . Whether or not this pos- 
sible cancelation would result in a net effective cosmo- 
logical constant even closer to the observational value is 
an intriguing issue for future investigation, where further 
interaction effects ought to be considered as well. 

We conclude with the following comment. We have 
shown that general relativity permits spontaneous fluc- 
tuations of metric, and the conformal fluctuations in par- 
ticular. These bear dynamical consequences on the evolu- 
tion of the averaged metric. Once the spontaneously fluc- 
tuating conformal field is set to regulate the net amount 
of vacuum energy, it appears to have somehow a life on 
its own. In this sense the conformal field is as real as the 
vacuum energy. One then wonders why we never have 
detected it. In this respect we notice that even standard 
gravitational waves are so weak that, to date, still fail to 
be detected. Hopefully projects such as LIGO and LISA 
will allow a direct unambiguous detection of GWs. In 
the case of vacuum conformal fluctuations, we are talk- 
ing about tiny modulations of space time that happen at 
the random scale, i.e. > 10 -35 m. These modulations 
average out at the classical scale. This makes them ide- 
ally very difficult to detect directly, even though there is 
a possibility that matter wave interferometry could pro- 
vide experimental evidence in measuring the related de- 
coherence on massive, non-relativistic quantum particles 
[H, [6^]. At this stage, it is clear however that, once the 
role of the conformal fluctuations in regularizing the vac- 
uum energy of the universe has been clarified, an indirect 
evidence of their reality could be, after all, the fact that 
observations seem to hint at a vacuum which is compat- 
ible with Lorentz invariance at the classical scale and a 
positive, small, cosmological constant. 

The inclusion of spontaneous conformal spacctime fluc- 
tuations as a fundamental ingredient of the quantum vac- 
uum seems to yield promising results. More efforts will 
be needed to clarify their deeper, quantum, meaning and 
cast further light upon the significance and the mutual re- 
lation between, quantum vacuum, cosmological constant, 
as well as the generation mechanism and the ultimate na- 
ture of mass. 
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APPENDIX A: ENERGY DENSITY AND 
PRESSURE OF MATTER FIELDS USING THE 
FREE FIELD APPROXIMATION 

Here we estimate the matter fields related quantities 
p*^ and pm in full physical units. From their definitions 
in (|62|) and (f72|) . we need evaluate the three quantities 
p^, pt and py. These can in principle be obtained with 
quantum field theory However this is not a trivial task. 
Within the Standard Model, vacuum energy density is 
estimated to be given by, at least, three main contribu- 
tions Q : (i) vacuum zero-point energy plus virtual parti- 
cles fluctuations, (ii) QCD gluon and quark condensates 
(iii) Higgs field. 

Through the following calculation we want to get a first 
approximation estimate, by neglecting fields interactions 
and by describing the free-field configuration as a collec- 
tion of decoupled harmonic oscillators of frequency 

uj k = Vc¥ + roV, (Al) 

where k is the norm of the spatial wave vector k and m 
is some effective mass of the field quanta under exam- 
ination. This is quite accurate for the EM field but it 
is likely to be quite a crude approximation for e.g. the 
QCD sector. By doing so we are neglecting higher or- 
der contributions to the vacuum energy density as well 
as the nonlinear and strong coupling effects of QCD and 
a detailed treatment of the Higgs fields. Nonetheless we 
expect to obtain a meaningful lower bound estimate to 
the vacuum energy. 

We consider an arbitrary matter field of effective mass 
m within the Standard Model. The field is comprised in 
$ or T depending on whether it is m = or m ^ 0. 
We calculate a lower bound to the associated vacuum 
energy density by estimating the overall energy density 
resulting from the summation of the zero point energies 
of all frequency modes, up to the random scale cutoff set 
by I = XL p. 

We model each independent field component as a scalar 
field 4* with mass m and associated Klein-Gordon stress- 
energy tensor 

T a K b G [4>] =4> a fo- \r, ab (W + ^0 2 ) ■ (A2) 

From this we see that (f> 2 has the dimension of force. Re- 
garding <j) as in its zero-point fluctuating state, we shall 



assume it to be stationary over spacetime and having 
an isotropic spectral density S^i^k), so that the mean 
squared field is given by 

^ = jirf f d3kS *W- (A3) 

Furthermore, the mean squared time derivative <j> t t = 
cififl satisfies 

<^>=(2^)3 Jd 3 kLo 2 k S4k) (A4) 

and the mean squared gradient (V0)i := <j) t i for i = 1, 2, 3 
satisfies 

<I V ^| 2 >=(^}3 / ' d'kk'S^k). (A5) 
It follows from (fAT]) . (fA3l) . dA4l) . (Mb that 

(Wo) = + <IW| 2 ) =-^f (|0| 2 ) . (A6) 

Using (|A4|) . (|A5|) and (|A6|) we see that the mean stress- 
energy tensor given by (|A2|) becomes 

(T^ G ) = (4>a<t> b ) , (A7) 
which yields the effective energy density 

P = (Too) = ^(4> 2 t ) = ^ysjd 3 k ^S^(k) (A8) 
and the effective pressure 

P= ^ Ti ') = l(l V ^ 2 > = ^)3 Jd 3 kk 2 S4k). (A9) 

A useful combination follows from (lATI) . dA8l) and dA9l) 

as 

p - ip = (2^)3^ / d'kS^k). (A10) 

The spectral density S^(k) itself can be determined 
through the well-known zero point energy density expres- 
sion 

"=(^?/ A ^ < A11 » 

that adds up contributions from all wave modes of (f>. 
Comparing (|A8|) and (|A11|) we see that the spectral den- 
sity S^k) takes the Lorentz invariant form 

he 2 

S*(k) = — . (A12) 

Substituting ([AT]) and (|Al"2j) into (fAUj) and ([510]) and 

integrating using J d 3 k = 4ir J dkk 2 up to the cutoff value 
k\ = jj 5 -, we obtain 
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and 



he 



,2 f k 



12tt 2 



dk k 4 



r> 2 



(A14) 



Introducing the dimensionless variable y := — , we can 
rewrite (|A13[) and (|A14|) as 



and 



12tt 2 V A J 



o v 71 + V 2 



(A16) 



in terms of the Planck energy density pp and the effective 
mass of the field in units of Mp/X, i.e. m\ := ( A j>/A) • 
For m\ <C 1 we can approximate (|A15|) . (|A16|) by 



2 / \ 2 

Pp Pp ( m 



and 



P 



A 4 4A 2 V M 



3A 4 12A 2 \M P J 



(A17) 



(A18) 



up to adding 0(m\) terms. This approximation is phys- 
ically well justified. The heaviest particles in the Stan- 
dard Model are the quark top, with m t « 173 GeV, the 
weak bosons , with myy± w 80 GeV and the Z boson 
with mz ~ 91 GeV. All the other particles have m < 1 
GeV. Being the Planck mass of the order of 10 19 GeV, it 
will be m\ < 1 as long as the cutoff parameter satisfies 
A < 10 15 . This is a safe upper bound, as the stochastic 
classical conformal fluctuations are expected to have a 
cutoff which should not exceed A tn 10 2 — 10 5 jUS i]. 
Within this approximation, by subtracting (|A17j) and 



„ Pp f m 



(A19) 



This relation can also be obtained directly from the right 
hand side of (|A10[) by following through the steps leading 
from (|AT3|) to (jAT8|) . 

The expressions (|A17[) and (|A18|) are the sums of two 
terms, one of which depends upon the effective mass of 
the field. Collecting the contribution from all massless 
fields we have the corresponding energy density c 2 p<a and 
pressure given by 



rrij and number of independent components Nx ■ Using 
(TAT7)) we have 



pr 



N r n 2 p P N 2 p P ( M V 



A 4 



4A 2 \Mp) ' 



(A21) 



where Ax := Y] a N-^j is the total number of independent 
components of the massive fields and 



M := 



(A22) 



is the weighted average squared effective mass. The total 
pressure of the massive matter fields then follows from 
(|AT8)) to be 



Pr 



NyTT 2 p P N r p P ( M 



3A 4 

We therefore have 



12A 2 \M P 



(A23) 



1, „ , NrppfM 



(A24) 



corresponding to the effective cosmological constant 



Am = 



NyAp ( M 



8A 2 V M 



(A25) 



Finally the leading contribution of p*^ can be obtained 
by substituting (|A20|) . (|A21)1 and ()A23|) into ([62]) to be 



, ir 2 Np P 

Pip ^4 ' 

up to adding 0(M 2 ) terms, where 



Ah 



M 



(Mp/X) ' 



and 



N := Nq, + N~t 



(A26) 



(A27) 



(A28) 



is the total number of independent components for mass- 
less and massive fields. The ratio between the traceless 
part and effective cosmological constant energy densities 
for the matter fields scales as 



e± oc mi « i, 

Pi, 



(A29) 



where the inequality holds within the Standard Model, 
with M w 10 2 GeV, and for any realistic value of the 
cutoff parameter. 



APPENDIX B: SOME TECHNICAL 
DERIVATIONS 



= 3p* = 



N^n 2 p P 
~Y 4 : 



(A20) 



where Nq, is the total number of independent components 
for all massless fields. 

In the massive case, the total energy density px follows 
from all contributing fields with individual effective mass 



The linearized Einstein tensor has the explicit form 



51 



GabVA ■=\d c dbK c + \d c d a h bc - \d c d c h ab 

- l -d a d b h - ^q ab (d c d d h cd - d c d c h) , (Bl) 
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where h := i] ab h ab denotes the trace. It is straightforward 
to verify that, for h ab = 2arj ab this yields 

G%> [2a Vab ] = 2i lab d c d c a - 2d a d b a, (B2) 

as used in (|55)) . 

Turning now to the calculation of (T^ b F ^) in (J49J) , the 

full explicit expression of — {G^j /8n can be simplified 

to yield the following nonlinear differential operator [5l[ 




(d a h cd d b h cd - x -d a hd h h - d a h bc d d h cd - dbKcddkA , 

(B3) 

where 

Kb ■= h ab - ^rj ab h (B4) 

is the trace reversed metric perturbation. Given this ex- 
pression it can easily be verified that, for h ab = 2r\ ab a, 

-^(^ ) [2H)--^(< 2 6 ) [H) = -^( W ). 

(B5) 



1(2) 

To conclude, the term E a j ) ' follows from taking the 
second order terms in equation ([7]). Expressing the co- 
variant derivatives using the linearized connection 

T ab ] ■= \v cd (daQM + d b q ad - d d q ab ) , (B6) 

where q ab '■= g a ]j is the first order effective metric per- 
turbation, gives 

8irZ ab 2) := (2V a a b -2g ab na) {2) 

= -rf' d (d a q bd + d b q ad - d d q ab ) a c 

+ VabV dCr l ef {d d q c f + d c q df - d f q dc ) a e 
~2q ab d c d c a + 2 Vab q dc d d a, c . (B7) 

Using q ab = /3 ab — 2ur\ ab to eliminate q ab , equation (|B7|) 
gives a series of terms involving products of f3 ab and 
a, whose corresponding averages vanish thanks to (|1 1|) . 
What is left is a sequence of terms, quadratic in a, which 
are obtained from (|B7|) upon substitution of q ab with 
— 2ar] ab . Their average can be evaluated to yield the 
final result 

8tt (^l b 2) ^=4 (a a a, b + 2 Vab a> c a c ) . (B8) 
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